UNCLASSIFIED 


AH  276  978 


luf  tUi 


ARMED  SERVICES  TECHNICAL  INFIMAniW  AGENCY 
ARLINGTON  HALL  STATIffll 
ARLINGTd  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE;  When  government  or  other  drawings^  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  irtiatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  ri^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


LZ-  3  -L 


00 


a 

<o 


ACIC  TECHNICAL  REPOKT  NO.  96 


PRINCIPLES  OF  ERROR  THEORY 
AND 

CARTOGRAPHIC  APPLICATIONS 


FEBRUARY  1962 


Aeronouticol  Chort  and  Information  C«ni«r 
United  Statoi  Air  Force 
St.  louit  18,  V.o. 


PRraCIPLES  OF  ERROR  THEORY 
ANL  C/J?rOORAHIIC  APPLICATTONS 

ACJC  Technical  Report  No.  96 
February  I962 


Prepared  by 
Clyde  R.  Greer.walt 
Melvin  E.  Shultz 


SUBI-IITTED: 


REVIE'/ED: 


THO-IAS  C.  FINN: 
Technical  Director 


APPROVED; 


IAJ' 


^  STEBBINS  W.  GRIFFITH  lont-l,  'JSAF 


Commander 


AERONAOTICAD  GHAHT  AITO  INFORMATION  CENTER 
Second  and  Areeiuil  Streets 
St  Louis  18,  Missouri 


NOTICES 


This  report  Is  issued  to  disseminate  timely,  useful 
technical  data  to  activities  engaged  in  cartography, 
geodesy,  and  related  subjects.  Nothing  in  this  report 
is  to  be  construed  as  necessarily  coinsidiug  with  USAF 
doctrine. 

This  report  does  not  contain  Iniormatlon  or  material  of 
a  copyrighted  nature,  nor  is  a  copyright  pending  on  any 
portion  thereof. 


ii 


CONTEMTS 


Page 


LIST  OF  ILLUSTRATIONS  v 

LIST  OF  TABLES  vl 

PREFACE  vil 

ABSTRACT  viii 

1.  ONE-DIMENSIONAL  (LKJEAR)  ERRORS  1 

1.1.  Introduction  1 

1.2.  Clasaes  of  Errors  1 

1.3.  Basic  Concepts  of  Probability  2 

1.4.  The  Normal  Distribution  of  a  Continuous  Random  Variable  6 

1.5.  Application  of  the  Probability  Density  Function 

to  Random  Errors  9 

1.6.  Precision  Indexes  l4 

1.7.  Conversion  Factors  17 

1.8.  Propagation  of  Errors  I8 

1.9*  Examples  of  Linear  Errors  2j 

2.  'RVO-DIMENSIONAL  (ELLIPTICAL,  CIRCULAR)  ERRORS  21 

2.1.  Introduction  21 

2.2.  Elliptical  Errors  23 

2.3.  Circular  Errors  23 

2.3.1-  Circular  Probability  Distribution  Function  23 

2.3.2.  Circular  Precision  Indexes  25 

2.3.3.  Circular  Conversion  Factors  4l 

2.3.4.  Propagution  of  Circular  Errors 


ili 


3.  THREE-DlMKiSrONAL  (ELLIPSOIDAL,  SPHERICAL)  ERRORS  it3 

3.1.  Introduction  U3 

3.2.  Ellipsoidal  Errors  1»3 

3.3-  Spherical  Ernrs 

3.3- 1-  Spherical  Probability  Distribution  Function 

3.3- 2.  Spherical  Precision  Indexes  4$ 

3.3.3-  Spherical  Conversion  Factors  52 

3.3.'+.  Propagation  of  Spherical  Errors  53 

4.  APPLICATIOM  OP  ERROR  THEORY  TO  POSITIONAL  INFORMATION  54 

4.1.  Positional  Errors  5^ 

4.2.  The  Accuracy  Statement  55 

4.3.  Summary  of  Formulas  and  Conversion  Factors  58 

APPENDIXES 

A  Percentage  Ihrobability  for  Standard  Error  Increments 
B  The  Moat  Probable  Value 
C  Propagation  of  Errors 

D  Derivation  and  Solution  of  the  Two-Dimensional  Probability 
Distribution  Flmctlou 

E  Derivation  of  the  Spherical  Probability  Distribution 
Function 

F  Substitution  of  the  Clrculai*  Form  for  Elliptical  Error 
Distributions 

REFERENCES 


iv 


ILLUSTRATIONS 


Figure  No.  TltJo  Page 

I  Probabilities  of  Ntunliers  from  the  Roll  of  Two  Dice  5 

2a  Normal  Probability  Density  Curve  8 

2b  Normal  Probability  Distribution  Curve  8 

3  Normal  Probability  Density  Curve  of  Ooserveo  Values  12 

It  Normal  Probability  Density  Curve  of  Errors  13 

5  Probability  Areas  l6 

6  Normal  Circular  Distribution  33 

7  Curve  of  the  P(R)  Function  When  P(R)  -  39- 35^  35 

8  Curve  of  the  P(R)  Function  When  P(R)  «  50.00?i  37 

9  MSPE  Probability  Curve  38 

10  Graph  of  Conversion  Factors  for  39»3^i>  to 

Probability  39 

II  Graph  of  Conversion  Factors  for  to  39*./^ 

Probability  Ito 

12  Comparison  of  Spherical  Standard  Error  Approxi¬ 
mation  Methods  50 

13  MRSE  Probability  Curve  51 


V 


TABLES 


Table  No.  Title  Page 

I  Suinmary  of  Linear  Precision  Indexes  IT 

II  Linear  Error  Conversion  Factors  17 

III  Values  of  the  Constant  K  23 

IV  Solution  of  P{R)  Function  for  P(R)  «  39  •3^'?’  3*^ 

V  Solution  of  P(R)  Function  for  P(R)  ■  ^0.00$  36 

VI  Summary  of  Circular  Precision  Indexes 

VII  Circular  Error  Conversion  Factors  ^1 

VIII  Values  for  the  Constant  W 

IX  Solution  of  P(S)  Function  for  P(S)  »  JO.OOjt  ^9 

X  Summary  of  Spherical  Precision  Indexes  52 

XI  Spherical  Error  Conversion  Factors  52 


Vi 


PREFACE 


Optimum  utilization  of  ACIC  research  and  production  requires 
that  the  accuracy  of  source  material,  interim  and  final  products  he 
considered.  The  accuracy  is  expressed  hy  an  error  statement  which 
indicates  whether  the  product  is  reliable  and  acceptable  or  should 
be  used  with  discretion.  Therefore,  the  error  statement  must  be 
representative  of  the  product  and  have  a  sound  statistical  basis. 

The  purpose  of  this  paper  is  to  present  and  explain  the  theory 
and  procedures  for  providing  a  valid  and  meaningful  orre-  statement. 

The  normal  distribution  of  linear  errors  is  explained  in  de¬ 
tail  because  two  and  three-dimensional  error  distributions  are 
more  easily  analyzed  statistically  ,  individual  treatment  of  the 
linear  components.  The  principles  of  the  linear  error  distribution 
apply  only  to  independent  random  errors,  assuming  that  systematic 
errors  have  been  eliminated  or  reduced  sufficiently  to  permit  treat¬ 
ment  as  random  errors. 

Although  a  truly  circular  or  spherical  error  distribution  seldom 
occurs  in  a  sample  of  observations,  the  concepts  are  desirable  for 
ease  of  compnutation  and  undersl-anding.  Consequently,  considerable 
atoe.'ition  is  given  to  the  computation  of  an  approximate  circular 
or  spherical  error  distribution  fr'ra  unequal  linear  componei.us  of 
a  two  or  three-dimensional  error  distribution,  yet  retaining  prop¬ 
erties  such  as  precision  indexes  of  the  truly  circular  and  spherical 
distributions.  Some  characteristics  of  circular  and  spherical 
error  distributions  differ  from  those  of  the  linear  error  distri¬ 
bution;  however,  the  distinction  is  of  an  academic  nature  and 
hence  is  not  sanphasized  in  the  text. 

Organizations  using  ACIC  charting  products  should  find  the 
discussion  helpful  in  interpreting  statements  of  cartographic  accu- 
I’acy.  The  formulas  and  principles  can  also  be  applied  to  weapon 
system  accuracy  evaluation  and  other  purposes  provided  that  the 
assumption  of  a  normal  distribution  of  independent  random  variables 
is  feasible. 

Important  functions  and  equations  are  presented  in  the  text, 
while  lengthy  derivations  are  relegated  to  appendixes.  Liberal 
numbers  of  references  are  inserted  after  major  headings  to  facilitate 
i'urther  study. 


ABSTRACT 


One  of  the  nost  useful  contributions  of  error  theory 
is  ‘che  precision  index  vhich  identifier,  the  error  distribution 
and  specifies  the  probability  that  the  true  error  in  a  quantity 
does  not  exceed  a  certain  value.  This  situation  is  applicable 
to  the  evalaauion  of  map  and  geodetic  information,  in  that  i^ 
maxes  possible  .meaniasful  accuracy  statements  having  uniform 
interpT'etation,  and  is  compatible  with  established  map  accui-acy 
standards  which  specify  limits  of  permissible  error  in  various 
categories.  Standardized  procedures  and  supporting  theory  iv-i 
computing  linear,  circular,  and  spherical  precision  indexes 
are  presented.  The  recommended  procedure  for  computing  the 
circular  or  spherical  standard  error  from  linear  standard 
errors  in  X  and  t,  or  X,  '£,  and  Z  directions,  respectively, 
is  to  average  the  linear  standard  errors.  Other  precision 
indexes  in  the  same  error  distribution  are  easily  computed 
from  the  linear,  circular,  and  spherical  stiindard  errors  — 
the  most  important  precision  indexes. 
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1.  ONE-DI^lt,’NSIONAL  (LINEAR)  ERRORS 


1.1.  Inbmluctlon.  Various  aspects  of  the  sciences  of  geofles.v, 
cartography,  and  photogranmietry  involve  the  measurement  of  physical 
quantities  and  the  utilization  of  such  measurements.  Regardless  of 
the  precision  of  the  instrument,  no  measurement  device  or  method 
gives  the  true  value  for  the  quantity  measured.  Mech/irj>.ai  impel ■■ 
factions  in  instruments  and  the  limitations  introduced  by  human 
factors  are  such  that  repeated  measurements  of  the  same  quantity 
result  in  different  values.  Variations  among  successive  values 

are  caused  by  errors  *  in  the  observations. 

While  the  theory  of  errors  does  not  yield  a  true  value 
nor  improve  the  quality  of  observations,  it  does  provide  a  vay  of 
estimating  the  most  probable  value  for  the  quantity  and  of  deter¬ 
mining  the  certainty  attributable  to  the  estimate.  Once  this  has 
been  established,  a  least  squares  adjustment  can  be  used  to  remove 
or  distribute  the  observational  errors  to  obtain  a  solution  which 
is  relatively  free  of  discrepancies. 

1.2.  Clasces  of  Errors,  (ref.  6,  19,  22)  Errors  fall  into 
three  general  classes  which  may  be  categorized  by  origin  as 

(l)  biuiiucrs,  (2)  systematic,  and  (3)  random. 


•The  true  error  of  each  observation  is  the  difference  between 
the  true  value  of  a  quantity  and  the  measui..(?.  value. 


Blvindex  s  are  mistakes  caused  by  misreading  scales,  trans¬ 


posing  figures,  erroneous  computations,  or  careless  observers.  They 
are  usually  large  and  easily  detected  by  repeated  measurements. 
Systematic  errors  follow  some  fixed  law  and  are  generally  constant 
in  magnitude  and/or  sign  within  a  series  of  observations.  T'he 
origin  of  systematic  errors  in  geodetic  measurements  is  primaiily 
within  the  instrument  or  measuring  device.  Causes  of  systematic 
error  Include  faulty  instrument  calibration,  errors  inherent  in  the 
graduation  of  scales,  and  changes  in  performance  resulting  from 
variations  in  temperature  and  humidity.  Systematic  errors  can  be 
eliminated  or  substantially  reduced  when  the  cause  is  known.  Random 
errors  are  tnose  remaining  after  blimders  and  systematic  errors  have 
been  removed.  Tney  result  from  accidental  and  unknown  combinations 
of  causes  beyond  the  control  of  the  observer.  Random  errors  are 
characterized  by:  (l)  variation  in  sign  —  positive  and  negative 
errors  occurring  with  equal  frequency,  (2)  small  errors  occurring 
more  frequently  than  large  errors,  and  (3)  extremely  large  errors 
rarely  occurring. 

The  probability  that  a  random  error  will  not  exceed  a 
cer+aln  magnitude  may  be  inferred  from  an  analysis  of  the  normal  or 
Gaussian  distribution  of  infinite  number  of  random  variables. 

1.3.  Basic  Concepts  of  Probability,  (ref.  2,  3)  probability 
is  defined  as  the  frequency  of  occurrence  in  proportion  to  the  num¬ 
ber  of  possible  occurrences,  or  simply,  the  ratio  of  the  niunber  of 
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successes  to  the  number  of  trials.  Let  A  and  B  symbolize  two 
completely  independent  events.  Denote  P(A)  as  the  probability  of 
the  event  "A"  and  P(B)  as  the  probability  of  the  event  "B".  The 
pT-nbabllity  of  any  event  happening  must  be  between  zero  and  one.' 
That  is,  zero  probability  means  that  the  particular  event  will 
never  take  place,  and  u  probability  of  one  means  that  t'v!  partiv:- 
ular  event  wil3  occur  each  trial.  For  example,  the  probability 
of  rolling  the  number  7  with  a  single  die  is  0.0  (an  Impossible 
event),  but  the  probability  of  roDling  a  number  from  and  in¬ 
cluding  1  through  6  is  l.G. 

Rule  1.  The  probability  of  event  A  is  equal  to  or  greater  than 
0  but  ec’-al  to  or  less  than  1. 

0  -  P(A)  -  1 

Rule  2.  The  probability  of  a  failure,  or  the  probability  of  an 
event  not  occurring,  is  1  minus  the  probability  that  it  will  occur. 

;  -  P(A)  =  failure  of  event  A 

Rule  3<  The  probability  of  either  of  two  events  A  or  B  occurring 
is  equal  to  the  sum  of  their  individual  probabilities. 

P(A  or  B)  =  P(A)  P(B) 


’  Probabilitj  is  also  denoted  by  a  percentage. 


An  example  is  the  prohabiilty  of  either  a  3  or  occurring  on  the 
single  roll  of  a  die: 

P(3  or  4)  =  1/6  +  1/6  -  1/3 

Rule  4.  The  probability  of  two  events  occurring  simultaneously  is 
equal  to  the  product  of  their  individual  probabiliLies. 

P(A  and  B)  =  P{A)  •  P(B). 

An  example  is  the  probability  of  both  A  «  3  and  B  »  4  occurring  in 
a  single  roll  of  2  dice: 

P(3  and  4)  =  1/6  •  i/6  =  1/36 
The  probabilities  of  occurrence  of  the  numbers  summed  from 
each  of  3^  possible  combinations  resulting  from  the  single  roll  of 
two  dice  are  presented  in  Figure  i.  The  probability  of  rolling  the 
number  7,  for  example,  is  6/36  or  I/6  since  there  are  six  combi¬ 
nations  which  have  a  sum  of  seven.  A  histogram'  of  the  data 
approximates  the  area  under  a  superimposed  smooth  curve.  If  the 
number  of  dice  in  a  single  roll  were  increased,  the  histogram 
would  rapidly  approach  the  smooth  curve,  called  the  normal  proba¬ 
bility  density  curve. 

*  A  column  is  constructed  for  each  number  by  blocks  each 
representing  an  area  equal  to  I/36  probability. 


Number 


Pi'Ouabillty  ' 


Combinations 


1  0 


2 

l/'36 

(1,1) 

3 

2/36 

(1,2)  {2,i) 

4 

3/36 

(1,3)  (3,1)  (2,2) 

5 

V36 

(1,4)  (4,1)  (2,3)  (3,2) 

r- 

5/36 

(1,5)  (5,1)  (2,4)  (4,2)  (3,3) 

1 

6/36 

(1,6)  (6,1)  (2,5)  (5,2)  (4,3)  (3,4) 

8 

5/36 

(2,6)  (6,2)  (3,5)  (*'.3)  ('*,'+) 

9 

Vi6 

(3,6)  (6,3)  (4,5)  (5,4) 

10 

3/36 

(i^,6)  (6,4)  (5,5) 

11 

2/36 

(5,6)  (6,5) 

12 

1/36 

(6,6) 

13 

0 

Probability 


6/36 


23456  769  10  11  12 


Ficure  1 

Probabilities  of  Numbers  from  the  Roll  of  Two  Dice 


*  Note  that  the  sum  of  the  probabilitj es  is  1. 
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l.*».  The  Normal  Distribution  of  a  Continuous  Random  Variable, 
(ref.  2*^)  The  area  under  the  normal  probability  density  curve 
(Figure  2a)  represents  the  total  probability  of  the  occurrence  of 
the  continuous  i-andom  variable  x  and  is  equal  to  one,  or  100%. 

The  mathematical  expression  of  the  curve  is  the  normal  probabi.lity 
density  function,  p(x): 

(x  -  n)" 

p(x)  =  e  20^  (i-l) 

o  -Jl* 

where:  x  =  the  random  variable 

4  =  a  parameter  representing  the  mean  value 
of  X 

a  =  a  parameter  representing  the  standard 
deviation,  a  measure  of  the  dispersion 
of  the  random  variable  from  the  mean, 

(I.  (The  square  of  the  standard  devia¬ 
tion  is  called  the  variance . ) 

'/is  -  2.5066  .  .  . 

e  •»  the  base  of  natural  logarithms,  2.71828... 
The  p.j’ameters  are  computed  from  an  infinite  number  of  random 
variables: 

4  -  — -  (1-2) 

n 


6 


a 


(1-3) 


i=l 


where : 

n  =  the  number  of  random  variables^  and 

n— *-  m  . 

The  normal  probability  distribution  function  determines 
the  probability  that  the  random  variable  will  assume  a  value 
within  a  certain  interval  and  is  derived  from  the  normal  proba¬ 
bility  density  function  by  integrating  between  limits  of  the 
desired  interval-  Letting  the  limits  range  from  -  •  to  x: 

.X 

p(x)  dx 


P(x)  =i 


aJTn 


(X  -  n)‘ 
20^ 


dx 


(1-4) 


The  value  of  Kx)  ranges  between  0  and  1,  illustrated  in  Figure  2b. 
As  x  approaches  its  upper  limit,  P(x)  approaches  1;  as  x  approaches 
its  lower  limit,  P(x)  approaches  zero.  This  is  true  since  x  cannot 
exceed  nor  be  less  than  its  defined  limits. 
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Figure  2a 

Nomal  Probability  Density  Curve 


P(x) 
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1.5*  Application  of  the  Probability  Density  Function  to  Random 
Errors,  (ref.  3,  15,  21,  22)  The  normal  probability  density  cvirve 
of  an  Infinite  number'  of  measurements  of  the  unknown  quantity  X  Is 
expressed  by  parameters  analogous  to  those  of  equation  (l-l).  The 
true  value  Is  the  metn  of  the  distribution  of  the  observed 
values  Xp  Xg,  X^  •••  Xj^.  The  curve.  Illustrated  In  Figure  3,  *.^3 
the  mathematical  form: 


(Xi  - 


p(X) 


2a 


a  \'2* 


(1-5) 


where : 


■V 


/r  (X,- 


Hv)' 


The  normal  probability  density  curve  of  errors  has  a  mean  of  zero 

and  is  Identical  in  form  to  that  of  the  observed  values.  Illustrated 

In  Figure  the  curve  Is  described  by  the  function: 

2 


p(€) 


*/ijr 


(1-6) 


where : 


the  true  error; 

€  -  -  Hx 


the  standard  deviation  of  the  errors,  here¬ 
after  designated  the  standard  error; 

r  .i® 

-2  _  ^1-1 


‘  The  population  or  universe  In  statistics. 
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Since  -chs  tn-e  value  of  a  quantity  cannot  "be  measured  and 
an  infinite  number  of  measurements  is  impractical,  estimated  values 


obtained  from  a  finite  number  or  sample  '  of  measurements  must  be 
substituted  for  the  true  value  and  the  parameters  of  the  density 
function.  The  most  probable  value  (X)  approximates  the  true 
value  and  is  determined  from  the  arithmetic  mean*  of  observed  values: 


z:.. 


(1-7) 


The  true  srror  in  approximated  by  the  residual  "x"*,  hereafter 
designated  the  error  and  defined  as  the  difference  between  the 
observed  value  and  the  most  probable  value: 


X  »  -  X 

The  standard  error  computed  from  a  sample  {a^)  is  identified  by 


(1-8) 


a  subscript* and  computea  from: 


4 


iV  (x^  - 

^i»l 


if 


n-1 


4E 


(1-9) 


•As  the  number  of  measurements  in  the  sample  becomes  larger,  the 
reliability  of  the  estimate  increases.  Often,  30  values  provide  an 
adequate  estimate. 

‘See  Appendix  B. 

‘The  residual  is  represented  by  "v"  in  some  texts. 

•  The  standard  error  derived  from  a  sample  is  designated  in  some 
texts  by  "s"  or  "m". 
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The  nomal  probahllity  density  function  of  errors  now  becomes; 


2 

X 

1  2a  2 

p(x)  =  - —  e  *  (l-io) 

°x 

The  parameters  X  and  may  assume  different  values  as 

various  samples  are  selected  from  the  same  population  and  are, 

therefore,  random  variables  with  dispersion  expressed  by  similar 

parameters.  The  standard  error  of  the  mean,  o_  >  and  the 

X 

standexd  error  of  the  standard  error,  ,  Indicate  the 

reliability  of  the  estimate  and  help  "round  off"  the  computed 
values :  , 


»  X  -  -  Ojj  X  +  Ojj  •  X  + 

Sample 


Figure  3 

Normal  Probatil^ty  Density  Curve  of  Observed  Values 
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-  € 

(Xi  <  4x) 


-  0  0  +  0 

(Xi  . 

Population 


+  € 

(Xi  >  Mjf) 


r 


p(x) 


(Xj^  <  X)  (Xj  -  X)  (Xi  >  X) 


Saunple 


Figure  4 

Normal  Probability  Density  Curve  of  Errors 
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1.6.  Precision  Indexes,  (ref.  3,  22)  A  precision  index  reveals 
how  errors  ai-e  dispersed  or  scattered  about  zero  and  reflects  the 
limiting  magnitude  of  error  for  various  probabilities.  For  example, 
505^  of  all  errors  in  a  series  of  measurements  do  not  exceed  ±20 
feet;  90^  do  not  exceed  ±49  feet.  Although  different  errors  are 
given,  each  expresses  the  asaie  precision  of  the  measuring  procer:' 
(Figure  5).  The  standard  error  and  average  error  (ij)  are  two 
Indexes  with  theoretical  derivations.  Common  usage  has  Included 
three  additional  probability  levels  which  are,  in  effect,  preci¬ 
sion  Indexes:  (l)  probable  error  (PE),  (2)  map  accuracy  standard 
(mas),  and  (3)  the  three  sigma  error  (3a). 

The  standard  error  is  the  most  Important  of  the  indexes 
and  has  the  probability  of; 


Hx)  » 


p(x)  dx  ■  0.6827 


(1-13) 


Or,  68. 27^  of  all  errors  will  occur  within  the  limits  of  ±  a^^. 

9Ttor  is  defined  as  the  mean  of  the  sum  of  the 
absolute  values  of  all  errors: 


V"  |(Xi  -  X)|  y 

^1-1  L  |X| 


n  n 

The  probability  represented  by  the  average  error  is  0.5751,  or 
57. 51^.  The  average  error  is  easily  computed  from  the  standard  error: 


0.797y  ffv 
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n 


(1-15) 


The  three  sigma  error,  as  the  name  implies,  is  an  error 
three  times  the  magnitude  of  the  standard  error.  The  3a  error  is 
used  because  it  approaches  near- certainty  —  0.9973  or  99*73?o 
probability. 
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1.7.  Conversion  Factors,  (ref.  20,  2?)  ^Ince  all  precis!  .i 
indexes  are  related  to  the  standard  error  (Table  l),  factors 
computed  from  this  relationship  (Table  II)  will  convert  the 
error  at  a  given  probability  to  the  error  at  another  probability. 


Table  I 

Summary  of  Linear  Precision  Indexes 


Symbol 

Probability 

Derivation 

PE 

.5000 

0.67^5  Ojj 

n 

.5751 

0.7979  oTx 

“^x 

.6827 

1.0000 

IttS 

.9000 

1.6449  ffx 

3o 

.9973 

3.0000  Ox 

1.3.  Propagation  of  Errcro.  (ref.  5>  29)  A  quantity  Is 
computed  from  tvo  measured  quantities  a  and  where  f(a,b)  denotes 
a  function  of  a  and  b.  The  error  Af  of  i.e  affected  by  the  errors 
In  both  a  and  b:  ^  and  db.  Assuming  a  and  b  are  Independent,  and 
the  errors  da,  db  are  randomly  distributed,  the  combined  error  AT 
can  be  computed  by  the  general  equation: 


where ; 


0^  ■  the  standard  error  of  f 


(1-20) 


Og,  0^  *  the  standard  errors  of  a  and  b 
df  df 

— ,  —  ■  partial  derivatives  of  f,  with  respect  t~ 
da  db 

a  and  b. 

Application  of  the  general  equation  to  specific  conditions  produces 
the  following  rules*  for  the  function  f(a,b): 


Rule  1.  Sum  and  Difference 

f  "  (a  +  b)  orf  ■  (a-b) 


Of. 


(1-21) 


‘Derivations  in  Appendix  C. 
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Rule  2.  Product  of  Factors 


a“b^ 


!£ 

f 

Rule  3»  Simple  Product  or  Quotient 

f  ■  a'b  or  f  ■  a/b 


(1-22) 


Indexes  other  than  the  standard  error  can  be  used  to 
propagate  errors.  For  example,  using  Rule  1: 


U-23) 


(PE).  -  '^(PE)a^  +  (PE)^^ 

"f  ■ 

(MAS)^  •'\J  {ma)^  +  (MAS)^^ 
and  1 3a) J  iSo)^  +  (3a)t^ 


However,  note  that  the  index  must  be  consistent  throughout  the 
formula.  That  Is; 
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nf  (MAS)j^^  +  {3<t)^2 

etc. 

1.9.  Examples  of  Linear  Errors.  The  foregoing  discussion  demonstrates 
the  use  of  the  normal  distribution  in  the  aned-ysls  of  random  errors. 

There  are  numerous  opportu'dties  for  the  occurrence  of  random  variables 
in  cartographic  und  geodetic  work.  For  example,  the  base  lines  and 
measured  angles,  observed  lengths  of  lines,  elevations,  etc.,  resulting 
from  geodetic  triangulation,  traverse,  and  leveling  all  contain  error. 

The  same  is  true  of  celestial  and  gravimetric  observations  as  well  as 
distances  measured  by  trilateration.  The  principles  of  error  theory 
can  be  used  advantageously  to  analyze  the  res'ilts  In  terms  of  the 
specifications  established  for  the  survey. 

In  ACIC,  the  normal  linear  error  distribution  has  Important  appli¬ 
cations  with  respect  to  evaluating  the  accuracy  of  positional  informa¬ 
tion.  In  addition  to  the  one-dimensional  errors  which  exist  in  such 
ppsitional  data  as  elevations  above  mean  sea  level,  the  linear  error 
components  of  two  and  three-dimensional  positions  can  be  analyzed  by 
applying  principles  of  the  normal  linear  error  distribution.  The  follow¬ 
ing  sections  contain  discussions  of  the  utility  of  the  linear  standard 
error  for  analyzing  two  and  three-dlnensloncil  distributi ore . 
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2.  TWO-DIMENSIONAL  (ELUPTICAL,  CIRCULAR)  ERRORS 


2.1.  Introduction.  A  two-dimensional  error  io  the  error  in  a 
quantity  defined  by  two  random  variables.  For  example,  consider 
the  true  geographic  position  of  a  point  referred  to  X  and  Y  axes. 
Each  observation  of  tne  X  and  Y  coordinates  will  contain  the  errors 
"x"  and  "y"-  When  assumed  random  and  Independent,  each  error  has  a 
probability  density  distribution  of; 


2 

X 


(Tj.  '/iit 


and: 


p(y) 


Applying  Rule  U  of  Section  1.3>;  the  two-dimensional  probability 
density  function  becomes: 


p(x>y)  • 


1 

2 


2x  dy 


2  2 

JL  + 

2  2 

\  '^x'"  V 


(2-1) 


Rearranging  terms: 


p(x,y)  Ojjffy  2k  -  e 


2  2 

a  2  g  2 
”x  “y 
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Therefore ; 


2  In  [p(x,y)  ®y  M  ^ 


(2-2) 


For  given  values  of  p(x,y),  the  left  side  of  equation  (2-2)  Is  a 
constant 


Then: 


2  2 
vSJLJL 
-  2  +  2 
®X  ®y 


(2-3) 


For  values  of  n(x,y)  from  0  to«»,  a  family  of  equal  prohabllity  density 


ellipses  are  formed  with  axes  K  and  K  ay. 
When  oy  ■  equation  (2-2)  becomes: 
~2a^  In  jp(x,y)  2xj  -  x®  +  y^ 


(2-4) 


For  a  given  value  of  p(x^y),  the  left  side  of  equation  (2-4)  is  a 
constant  which  is  the  square  of  the  radius  of  an  equal  probability 
density  circle. 

The  probability  density  function  Integrated  over  a  certain  region 
becomes  the  probability  distribution  function  which  yields  the  proba¬ 
bility  that  X  and  y  will  occur  simultaneously  within  that  region,  or; 


(x,y)  -  I J v{x>y) 


dx  dy 


However,  since  both  positive  and  negative  values  of  either  "x"  or  "y" 
will  occur  with  equal  frequency,  the  errors  may  be  considered  as  radial 


errors,  designated  by  "r",  where  r  *  v  x2  +  y- 


2.2.  Elllptloal  Errorg.  (ref.  15^  20)  lije  probability  of  an 
ellipse  Is  given  by  the  distribution  function; 

K® 

P(x,y)  -  1  -  e"  (2-5) 

Hie  solution  of  equation  (2-5)  with  values  of  K  for  different  proba¬ 
bilities  yields  the  results  shown  in  Table  III.  For  a  395t  probability, 
the  axes  of  the  ellipse  are  1.0000  Ojj  and  1.0000  Oyj  for  a  ^0$  proba¬ 
bility,  the  axes  are  1.177^  Ox  and  1.177^  Oy. 


Table  III 

Values  of  the  Constant  K 


Probability 

K 

39.35i 

1.0000 

50.005^ 

1.177*^ 

63.21S( 

1.1^142 

90.00* 

2.i460 

99.00^ 

3.0349 

99.78?^ 

3.5000 

The  use  of  the  error  ellipse  Is  complicated  by  the  problem 
of  axes  orientation  and  propagation  of  elliptical  errors.  Therefore, 
the  ellipse  is  commoniy  replaced  by  a  circular  form  which  Is  easier 
to  use  and  understand. 

2.3.  Circular  Errors. 

2.3.1.  Circular  Probability  Distribution  Function,  (ref.  i, 
2k)  The  probability  distribution  function  of  the  radial  error 
expressing  the  probability  that  "r"  will  be  equal  to  or  less  than 
radius  R,  or  the  probability  that  the  vector  xy  vlll  be  contained 
within  a  circle  of  radius  R,  Is  derived  in  Appendix  D  and  stated  as: 
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[‘•a 


dr  (2-6) 


A  special  ease  of  the  P(P.)  function  (2-6)  Is  formed  vhen 


and  Appendix  D>  part  2: 


P(R)  -  Pg  -  1  -  e 


(2-7) 


where : 


P  •  the  circular  prohablllty  distribution  f<mctlon>  a  special 
c 

case  of  P(R) 

R  ■  the  radius  of  the  probability  circle 
Oq  "  the  circular  standard  error >  a  special  case  of  vhen 

•  «y 

When  Ojj  and  Cy  are  not  equal,  the  P(R)  function,  (2-6),  is 

Oy 

modified  by  letting  "a"  equal  the  ratio  —  where  is  the  smaller 


standard  error  of  the  two.  Then  from  Appendix  D,  part  3; 


A 

■  —f 


-V  ^  \  , 

e  VVA/  dv 


(2-8) 


‘  ‘  ^  [“J 


2k 


Equation  (2-6)  can  be  solved*  for  different  probabilities  or  values 
of  P(R)  representing  precision  indexes  of  the  error  distribution. 

2.3.2.  Circular  Precision  Indexes,  (ref.  19>  20,  27)  The 
precision  Indexes  Illustrated  in  Figure  6  are  measures  of  the  disper¬ 
sion  of  errors  in  a  distribution  and  represent  the  error  which  is 
unlikely  to  be  exceeded  for  a  given  probability.  The  preferred  cir¬ 
cular  precision  indexes,  consistent  with  indexes  used  in  the  linear 
distribution,  are;  (l)  the  circular  standard  error  (Og),  (2)  the 
circular  probable  error  (CPE,  CEP),  (3)  the  circular  map  accuracy 
standard  (CKAS),  and  (l^)  the  circular  near-certainty  error,  three-fl-'- 
sigma  (3.5  .  The  mean  square  positional  error  (MDPE),  on  additional 

index  which  has  been  used  at  ACIC,  is  not  recommended  because  the  proba¬ 
bility  represented  varies  when  and  Py  are  not  equal. 

The  probability  of  the  circular  standard  error  is 

found  by  solving  equation  (2-7)  for  Pg  when  Og  ■  R,  thus: 

2 


» Described  in  Appendix  D,  part  ^ . 
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_  1 

a  1  -  e  2 

-  1  -  0.60653 

..Pc  -  0.3935  (2-9) 

That  lE,  39. 35^  of  all  errors  In  a  circular  distribution  are  not 
expected  to  exceed  the  circular  standard  error. 

For  a  truly  circular  distribution#  the  linear 
standard  errors  are  equal  and  Identical  to  the  circular  stemdard 
error  (o^  "  ■  ®c^*  When  Ojj  and  Oy  are  not  equal#  a  normal 

circular  error  distribu.tion  may  be  substituted  for  the  elliptical 
distributior..  The  substitution  is  satisfactory  for  error  analyslp 
within  specified  ®^n/®max  *  ratios.  Because  of  distortion  in  the 
error  distribution  ’  for  low  ratios#  however#  the  circular  concept 
should  be  used  with  dis^reticn. 

An  approximate  circular  standard  error  is 
determined  from  equation  (2-8)  by  letting  P(r)  ■  39.355^  and  R  ■  Oq. 
Values  of  g^/Onmy  for  ratios  of  gniin/^max 

tained  in  Table  IV  and  plotted  in  Figure  7.  For  the  Oniin/'^max 
between  1.0  and  0.6#  the  curve  is  a  straight  line  with  the  equation; 


•  Where  minimum  or  smaller  linear  standard  error  of  the 

two. 

*  See  Appendix  F. 
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(0.5222  Onjin  Omax^ 


(2-10) 


a 


A  rapid  approximation  gives  a  slightly  larger  Oq  value  for  the  same 
®niin/®mex 

^  0.5000  (ojj  +  Oy)  (2-11) 

‘^min/*^max  approaches  zero,  the  39- 35?^  probability  curve  follows  a 
transition  from  circular,  through  elliptical,  to  the  linear  distri¬ 
bution  form.*  The  curve  does  not  effectively  represent  a  circuit 
standard  error  for  Otnln/®max  *“a^ios  less  than  u.6  because  it  is  not 
compatible  with  other  circular  precision  indexes.  For  example,  the 
factor  1.177^  convesTvS  a  circular  error  at  39?^  probability  to  a 
circular  error  at  30%  probability  when 

®mln  "  factor  converting  a  linear  error  at  39^  probability 

to  a  linear  error  at  30%  probability  is  1.309^.  *  The  circular  stand¬ 
ard  error  computed  from  equation  (2-11),  however,  can  be  converted 
to  other  circular  precision  Indexes  by  constant  circular  conversion 
factors  *  for  Ond.n/%ax  ratios  between  1.0  and  0.2  and  is,  therefore, 
the  preferred  method  for  approximating  the  circular  standard  error. 


•When  ■  0,  the  factor  0.5151  converts  a  linear  error  at  6856 
probability  to  an  error  at  39* 35^^  probability. 

•  The  transition  curves  of  conversion  factors  are  shown  in  Figures 
10  and  11. 

•Presented  in  Section  2. 3. 3. 
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»  0.69315  (2a/) 

R  «  1.177^^  a^j 

CPE  «  1.177^*  a  {2-13) 

c 
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When  Ox  and  Oy  are  not  eq.ual>  on  approximate  CFli.  la  determined  from 
equation  (2-8)  by  letting  P(R)  ■  JO.OOJJ  and  R  ■  CPE.  Values  of  CPE/cj^ 
for  ratios  of  «^n/®ma.t  tabulated  in  Table  V.  The 

50^(  probability  curve  plotted  In  Figure  8  is  approximated  by  a  series  of 
straight  lines  for  different  ratios  cf  djoin/^^max  equations: 

CPE  «  (0.6i42  +  0.5632  ffnax)  (2-l4) 

when  ®iBln/®max  between  1.0  and  0.3 

CPE  -  (0.4263  +  0.6196  0^)  (2-15) 

when  o  ._/o  is  between  0.3  and  0.2 
min  max 

A  rapid  approximation  of  the  CPE  plots  as  a  straight  line  wUch  inter¬ 
sects  the  50^(  probability  curve  at  the  point  where  *0.2  and 

has  the  equation: 


CPE  0.5887  (dx  +  Oy)  (2-l6) 

when  Oain/^max  between  1.0  and  0.2 

The  CPE  computed  by  equation  (2-l6)  is  compatible  with  the  circular 
standard  error  computed  by  equation  (2-11)  '  and  iS/  therefore,  the 
preferred  method  for  approximating  the  circular  probable  error  within 
the  specified  limits. 


‘  That  Is,  the  conversion  factor  of  1.1774  for  converting  Og  to 
CPE  is  conbtant  for  ratios  of  ajjjin/<%jax  between  1.0  and  0.2.  Nov.e  that 

1.1774  [0.5000  (ox  +  dy)1  -  0.5887  (Ox  +  ‘^y)* 
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Although  a  circular  error  concept  Is  not  reconmended  for  aoin/^mejc 
ratios  less  than  0.2,  a  i.ear-linear  JOjt  prohabllity  error  may  be 
computed  to  represent  a  CPE  for  lower  ratios  when  a  comparison  of 
circular  errors  derived  from  different  sources  is  required: 


CPE  ~ 

(0.2l4i  0.6621  Oojax) 

(^»-37) 

when  ffnin/Ojiax  between  0.2  and  0.1 

CPE  ~ 

(0.0900  +  0.6745  0^) 

(2-18) 

when  Ojjin/Omax  between  0.1  and  0.0 

CPE  ~ 

0.6745  o„«x 

(2-19) 

when  ° 


The  following  alternate  methods  of  computing  an  approximate  CPE  are 
not  recommended  because  of  limited  applicability: 


(2-20) 


and  CPE  ~  O.8325  +  0^  (2-2l) 

when  Oniin/^max  between  1.0  and  0.8 

The  mean  square  positional  error  (ref.  1,  11 )  is 
defined  as  the  radius  of  the  error  circle  equal  to  1.41^2  and  has 
little  significance  in  a  truly  circular  error  distribution.  However, 
when  Ox  and  Oy  are  approxlmaxely  equal,  the  MSP^J  defines  the  error  in 
a  gcogmpmc  position  and  io  computed: 
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=  63.215& 


(2-24) 


When  a„  4  the  solution  of  (2-23)  yields  values  of  P«  (plotted  in 
Figure  9)  ranging  from  64^  when  ajn^n/Oj^  =  0.8  to  77^  when  cfn^n/‘^max 
••0.3.  Because  of  the  variation  in  prohahllity,  the  MSPE  is  n't 
recommended  for  use  as  a  precision  index. 

The  circular  map  accuracy  standard  is  based  on  the 
percentage  level  in  use  by  the  U. S.  National  Mtip  Accuracy  Standards 
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vhlch  specify  that  no  more  than  105^  of  the  well-defined  points  in  a 
map  will  exceed  a  given  error.  The  standards  are  commonly  Interpreted 
as  limiting  the  size  of  error  which  90$  of  the  well-defined  points 
will  not  exceed.  Therefore,  the  circular  map  accuracy  standard  is 
represented  by  the  value  of  R  in  equation  (2-7)  when  P^,  ■  O.9O,  and 
is  computed: 


CMAS  -■*  2.1U60  Og 

(2-25) 

or 

CMAS  =  1.8227  CPE 

(2-26) 

The  three-five  sigma  error,  representing  a 
circular  probability  of  99*785^,  approaches  near-certainty  in  a  circular 
distribution  and  nas  a  magnitude  3.5  times  that  of  the  circular  st.*':’*- 
ard  error. 
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Table  IV 


Solution  of  P(R)  Function  for  P(R)  -  39 >35^ 


®inln 

®max 

1 

®mBX 

1.0000 

1.0000 

0.8165 

0.9063 

0.651*7 

0.8197 

0.5000 

0.7323 

0.3333 

0.6327 

0.2294 

0.5727 

0.1005 

0.5274 

0.0 

0.5151 

Note:  When  P(R)  •  39-35^/  R 


3I+ 


35 


0.50 


Table  V 


Solution  or  P(R)  Function  for  P(R)  -  50.00?t 


f 

I 
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2. 3*3*  Clr<‘.ular  Conversion  Fac:tors.  (ref.  20,  27)  The 
relationships  of  the  circular  standard  error  to  other  circular  pre¬ 
cision  indexes  are  summarl 7.ed  in  Table  VI.  Conversion  factors  (Table 
VII )  computed  from  these  relationships  convert  a  circular  error  ut  a 
given  probability  to  a  circular  error  at  another  probability.  When  a 
circular  error  dlstriijuticn  is  subSwituted  for  un  ellip'  U:al  ulwirllu- 
tion,  the  circular  conversion  factors  are  retained. 


Table  VI 

Summary  of  Circular  Precision  Indexes 


Symbol 

Probability 

Derivation 

®c 

.3935 

1.0000  Og 

CPK,  CEP 

.5000 

MSPE 

.6321 

1.4i1^2  Oj, 

CMAS 

.9000 

2.]46o  Og 

3.5  ffc 

.9978 

3.5000  Og 

Table  VII 


Cirrular  Error  Conversion  Factors 


Froni-''..^ 

39.35% 

30.00% 

63. 21^ 

90.00% 

99.7Q% 

39.35^ 

30.ooi> 

63.21^ 

S0.00<f, 

99- 78^6 

1.0000 

O.8U93 

0.7071 

0.4660 

0.2857 

1.1774 

1.0000 

0.8325 

0.5486 

0.3364 

1.4i42 

1.2011 

1.0000 

0.6590 

0.4o4o 

2.l46o 

1.3227 

1.5174 

1.0000 

0.6131 

3.5000 

2.9726 

2.4749 

1.6309 

1.0000 

I-.  Propag.ation  ol*  Circular  Errors,  (ref.  5,  29)  A  two- 
dimensional  quantity  derived  from  a  n’jmbcr  of  independent  variables  has 
a  circular  error  resulting  from  tne  errors  in  each  variable.  The  total 
circular  error  is  determined  by  propagating  the  linear  components  in 
each  of  tne  two  dimensions  by  methods  described  in  Section  1.8.,  and 
computing  the  circular  form  by  the  procedure  shOTm  in  Section  2.:!. 2. 

For  example,  the  total  circular  error  of  a  quantity  C^,  derived  from 
C,p  ■  Cl  f  Cg  +  ...  Cj^,  is  found  by: 


a 


% 


0.5000  (a„  +  ) 


(2-27) 


An  alternate  approximate  propagation  method  combines  the  circular  error 
of  each  Independent  variable  directly,  thus; 


(2-28) 


Precision  indexes  other  than  the  standard  error  may 
be  used;  however,  the  index  must  be  consistent  throughout  the  computs- 


tlons. 


3.  THREE-DIi.lENSIOML  (ELLIPSOIDAT,,  SPHERICAL)  ERRORS 


3.1.  Introduction.  A  three-dimensional  eri’or  is  the  error  in  a 
quantity  defined  by  three  random  variables.  Expanding  on  the  example 
in  Section  2,2.,  a  point  is  referred  to  X,  Y,  and  Z  axes  which  estab¬ 
lish  the  spatial  position  of  the  poi.nt.  Wlien  random  and  independent, 
the  errors  x,  v,  and  a  “ac!'.  have  a  IMiear  probability  Jl.r.ributlon. 
The  three-dimensional  probability  density  function  is  expressed  by: 


p(x,y,z) 


.2 

Oy  0^ 


.  2  La  2 


JC 


(3-1) 


Similar  to  Section  2.1.,  the  probability  density  function  can  be  written; 


W 


(3-2) 


where :  ^ 

w2  =  -  2  In  ^p(x,y,z)  o^.  (2«)^j 

p 

For  values  cf  the  constant  W  from  O-*-  ob,  the  density  function  defines 
a  family  of  ellipsoids  of  equal  probability  density. 

3.2.  Ellipsoidal  Errors,  (ref.  15,  20)  The  probability  of  an  error 
ellipsoid  is  given  by  the  probability  distribution  function: 

-  -  t^ 

t^e  2  at  (3-3) 
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where : 


s  =  the  radial  error; 


+  V- 


Oj,g  =  standard  error  of  the  radial  error  "s" 

The  solution  of  equation  (3-3)  for  W  yields  the  values  given  in  Ta';le 

vm. 


Table  VITI 

Values  for  the  Constant  W 


Probability 

W 

19.9$ 

1.000 

50 

1.538 

6o.C 

1.732 

90 

2.500 

99 

3.368 

99.89 

4.000 

3.3*  Spherical  Errors. 

3 • 3 • 1 •  Spherical  Probability  Distribution  Function .  ( ref . 20 ) 
When  Ox  “  '^y  "  ‘^z  “  ®rs  ®  (3-l)  becomes  the  spherical  prob¬ 
ability  density  function: 


p(8) 


(3-M 


U4 


wl!i}''c :  0,^  ~  spl  erica]  stanc^ard  error 

Ir (.ograi-in;^  p(£)  from  s  =  0  t.o  3  =  S,  cojial  ion  (3"‘^)  Jicroincs  tiip  sp'iari- 
cu]  probal'ility  distribuLion  I'lnciion:" 


w)\ere:  S  -  radius  of  fue  probal  ill'.^/  spliere 

Equation  (3-;)  can  be  sojved  by  an  appi-oxination  fonrula  (ref.,  Ij, 


-13): 


P(S) 


1.253  -  Ce 


2 


C 


n 


(3-^) 


wbero C 


S_ 

0 

s 


3.3*2.  c.p'nc-rical  Precision  Indexes,  (ref.  20,  27)  A  sahe"]- 
cal  error  distrib.tion  is  represented  by  inde.xes  si'miJar  to  tiiose  in  Sec¬ 
tions  ].6.  and  2.3.2.  Prefoned  sp':cri<-al  recision  indexes  include: 

(1)  itie  spherical  standard  error  {o^),  (2)  f'  e  spherical  probai  io  error 
(SPE),  (3)  t'>e  spherical  accuracy  standard  (SAS),  and  (^)  the  snheri- 
cal  near-certainty  error,  four  sitir.a  (^Og).  The  -nean  radiaJ  spherical 
error  (I'O'iSE),  an  index  vfiich  has  been  used  at  AOTr’,  is  not  rccorir-ended 
because  ttie  uro'oabilLty  represented  vt.'-ies  wiion  a^,  0  ,  and  a  am  not 

X  y  z 

equal. 


See  Appendix  E. 


The  probability  of  an  error  splicrt;  oi  radius  equal 


to  the  spherical  standard  error  is  computed  by  equation  (3-6)  for  the 


condition  C  = 


1  as  follows;. 


0.7970846 

0.60653 


0.67032 


C.8e‘°’*^  =  0.53626 

P(S)  -  0.79788  (1.253  -  0,6065  -  0,3948) 

P(S)  0,20  or  20^  '  (3-7) 

For  a  truly  spherical  distribution,  the  linear  stand¬ 
ard  errors  are  equal  and  identical  to  the  spherical  standard  error 
(ff^  =  Oy  =  Oj  3  Og).  When  a^,  Oy,  and  are  not  equal,  the  spherical 

standard  error  is  approximated  by: 

05  ~  0,3333  (Ojj  +  Oy  +  a^)  (3-8) 

when  between  1.0  and  0.35 

The  substitution  of  a  spnerical  form  for  an  ellipsoidal  distribution  is 

not  recommended  when  the  0,„;j^„/*^n!ax  than  0.35' 

The  following  alternate  metiiod  of  approximating  Og 

Is  net  recommended  because  of  jimited  applicability:* 


'  A  more  accurate  value  is  determined  by  an  expansion  !n  series  to 
be  19 •9';®  probability. 


*  Figure  12  compares  curves  computed  from  equations  (3-8)  and  (3-9), 


(3-9) 


2  P  P 

V  ^  V  ^ 


when  o^in/'^max  between  1.0  anO  0.9 
The  ."pherjeal  probable  error  is  defined  as  ihe 
magnitude  of  the  spherical  radius  S  when  the  function  P{S)  =*0.5  or 
50^^.  Expressed  in  the  form  S  ■  C  a^,  the  spherical  ^.r' beVle  error  ib 
computed  by: 


SPE 


1.53B2  o 


(3-10) 


The  P{R)  function  for  two-dimensional  errors  is  solved  by  the  use  of 
Grad  and  Solomon's  tables.*  Expanding  this  method  into  the  spherical 
distribution,  the  radius  S  for  e  505^  probability  sphere 
computed  in  terms  of  for  ratios  of  amin/®max  ^mid/'^max 
tabulated  in  Table  IX.*  Utilizing  these  values,  an  approximation  o^ 
the  spherical  probable  error  can  be  computed:'  * 

SPE 0.5127  (ox  +  Oy  +  (3-^) 

when  Oinin/®max  between  1.0  and  0.35 

The  mean  radial  spherical  error  is  the  radius  of 
the  error  sphere  equal  to  1.732  Og,  or ‘IYoq,  in  a  truly  spherical 
distribution.  When  0^  ¥  <^y  ¥  <^z>  bhe  MRSE  is  computed  by: 


*  See  Appendix  D. 

’where:  =*  the  minimum  sigma,  or  smallest  standard  error  of 

the  tiiree, 

*^max  “  maximum  sigma,  and 
'^mid  middle  sigma. 

’Note  that  l.pjoii  [0.3333  (®x  ^  ”  0.5127  {a^  *■  Oy  o^). 


kl 


KRSE 


~  V? 


2  2 
,  ■*  a. 


(3-12) 


when  is  between  1.0  anb  0.0 

min'  max 


The  probabilities  representerl  by  the  MRSE  are  computed  by  equation 
(3-6).*  Because  of  the  vai-iation  in  probability,*  the  MRSE  is  not 
recommended  for  use  as  a  precision  index. 

The  spherical  accuracy  standard  is  defined  as 
the  magnitude  of  the  spherical  radius  S  when  the  function  P(S)  ■ 

0.9  or  90^.  Expressed  in  the  form  S  =  C  Og,  the  spherical  accuracy 
standard  is  computed  by: 

SAS  -  2.500  Og  (3-13) 


The  four  sigma  error,  representing  a  spherical 
probability  of  approaches  near-certainty  in  a  spherical 

distribution  and  has  a  magnitude  four  times  that  of  the  spherical 
standard  error. 


‘  Illustrated  in  Figure  13. 

•When  '1^®  probability  is  6o.825tj  when  =  10, 

Oy  *  3^  snd  a  6,  the  probability  is  69.36^f’. 
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Table  IX 


Solution  of.  F(S)  Function  for  P(s)  «  50.005t 


®raid 

°max 

®inln 

®max 

SPE-.  0.5127  (Ojj  +  ffy  +  cj 

%BX  “  ^ 

0.866 

0.866 

1.4007 

1.0 

0.707 

1.3892 

1.3879 

0.775 

0.632 

1.2341 

1.2341 

0.577 

0.577 

“.«0< 

1.1044 

o.agU 

1.2002 

0.707 

0.4o8 

1-0894  0^^ 

1.0844 

0.535 

0.378 

0.9791 

0.9808 

0.35'+ 

-j 

0.354 

_ 

0.8689 

0.8757 

Figure  12 


f 
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(10,10,10)  (10,1 
(10,7,9) 


3. 3- 3-  Sp^c.lcal  Conversion  factors,  (ref.  kO,  27)  The 
relationships  of  the  spherical  standard  error  to  other  spherical  pre¬ 
cision  indexes  are  summarized  in  Table  X.  Conversion  factors  (Table 
XI )  computed  from  these  relationships  convert  a  spherical  error  at  a 
given  probability  to  a  spherical  error  at  another  probability. 


Table  X 

Summai-y  of  Spherical  Precision  Indexes 


Symbol 

Probability 

Derivation 

®s 

.199 

1.000  Sg 

SPE 

.50 

1.539  Os 

MRSE 

.608 

1.732  Og 

SAS 

.90 

2.300  Og 

.99^9 

4.000 

9 

Tabic  XI 

Spherical  Error  Conversion  Factors 


FroiTr-^ — ^ 

19.956 

5056 

60.856 

9056 

99. 89^ 

19.95& 

1.000 

1.538 

1.73? 

2 . 500 

4,000 

50^ 

0.650 

1.000 

1.326 

1.625 

2.600 

60.856 

0.577 

0.888 

1.000 

1.443 

0 

CVJ 

90^ 

0.400 

0.615 

0.693 

1.000 

i.6co 

99.8956 

0.250 

0.385 

0.433 

0.625 

_ 1 

1.000 
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three-dimensional  quantity  derived  from  a  numhfr  of  independent 
variables  has  a  spherical  error  resulting  from  the  errors  in  each 
variable.  The  total  spherical  error  is  determined  by  propagating 
the  linear  components  in  each  of  the  tlii-ee  dimensions  by  methods 
described  in  Section  l.d.>  end  computing  the  spherical  foi.ti  by  the 
procedure  shown  in  Section  S-S-S*  For  example,  the  total  spherical 
error  of  a  quantity  S^,  derived  from  S,p  ■  +  Sb  +  . .  .  found 


0.3333  io..  t  a,  +  a,  ) 


(5-1^) 


An  alternate  approximate  propagation  method  combines  the  spherical 
error  of  each  independent  variable  directly,  thus: 


Precision  indexes  other  then  the  standard  error  may 
be  usedj  however,  the  index  must  be  consistent  th/oughout  the  co:iiputa- 


tions. 


4.  APPLICATION  OF  ERRO'^  THEORY  TO  POSITIONAL  INFORMA.TION 

4.1.  Positional  Errors.  By  the  use  of  error  theory  in  the  eval¬ 
uation  of  ACIC  positional  information,  it  is  possitle  to  establish  a 
meaningful  accuracy  statement  subject  to  .ini form  interpretation.  To 
provide  a  logical  and  acceptable  basis  for  computation  and  comparison, 
positional  errors  are  assumed  to  follow  a  normal  disti ibution.  Ti.t 
assumption  is  valid  because  positional  error  components  generally 
follow  a  normal  distribution  pattern  ^en  sufficient  data  is  available. 

The  statistical  treatment  of  errors  is  applied  to  measurable  quanti¬ 
ties  found  in  the  sources  of  positioning  information.  The  d'  fferences 
between  the  surveyed  coordinates  of  ground  control  and  the  scaled  co¬ 
ordinates  of  the  same  control  symbolized  on  maps  are  considered  to  be 
the  errors  in  the  geodetic  base  of  the  map.  Analysis  of  the  lineeir 
components  —  latitude  and  longitude  or  grid  Northing  and  Easting  — 
provides  a  two-dimensional  expression  of  the  accaracy  of  the  geodetic 
base,  'ftien  all  the  linear  standard  errors  occurring  dui-ing  map  con¬ 
struction  are  combined  and  converted  to  a  circular  distribution,  the 
final  map  accuracy  statement  is  expressed  in  terms  of  circular  errors. 

Among  the  positioning  errors  in  maps,  there  axe  often  those  which 
are  not  measurable  and  \diich  must  be  estimated  by  empirical  methods. 
When  this  is  necessary,  an  additional  assumption  must  be  made  to  the 
effect  that  such  data  is  compatible  with  computed  data  and  that  empiri¬ 
cally  derived  error  data  will  also  follow  the  tneorei>ictil  error  dis- 
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VaT.'ou'-  types  oT  points  require  different  p'U’amcters  to  establibh 
p-'-'rise  positions.  These  have  been  discussed  as  one,  two,  or  three- 
Jimenstona]  coordinates.  For  example,  a  vertical  position  (elevation) 
requires  only  a  one-diiJionjixnel  '•nordinate  —  the  height  of  t)ie  point 
above  a  reference  datumj  a  geodetic  position  Is  expressed  by  two- 
dimensional  coordinates  —  latitude  and  longitude  rcl tvi  ner-’  ir  « 
specific  datum;  and  spatial  positions  require  three-dimensional  coor¬ 
dinates  such  as  the  x,  y,  z  coordinates  in  a  rectangular  system.  "Tic 
errors  accumuiated  in  the  pri'ccss  of  determining  the  various  positions 
must  be  evaluated  in  the  same  dimensions  required  to  express  the  posi¬ 
tion.  Errors  for  vertical  positioning  can  be  assumed  to  follow  a  normal 
linear  distribution;  those  for  a  geodetic  position  —  a  circular  dis¬ 
tribution;  and  the  errors  for  a  spatial  point  can  be  assumed  to  follow 
a  normal  spherical  distribution. 

4.2.  'fhe  Accuracy  Gtatement.  Two  major  groups  of  data  fall  within 
Air  Force  positioning  requirements:  (l)  maps,  charts,  and  other  graphics 
and  (2)  specific  points.  By  the  use  of  error  theory,  a  horizontal  ac¬ 
curacy  evaluation  of  the  graphic  as  a  whole  can  be  obtained,  i.e.,  a 
specified  probability  that  the  true  errors  in  well-defined  planimetry 
will  noc  exceed  the  given  quantity.  Fiap  accuracy  can  also  be  inter¬ 
preted  as  percentage  —  the  percentage  of  well-defined  points  which 
will  not  contain  errors  exceeding  the  given  magnitude.  Similarly, 
vertical  accuracy  is  stated  as  a  given  probability  that  the  Hnear 
errors  in  vertical  nj.sitiou  tu’e  nxt  likely  to  eAceeu  a  sciecified  .’alue. 


The  accuracy  of  a  specific  point,  is  expressed  also  by  a  statement,  of 
probability  and  error  magnitude.  TTie  occiu'acy  statement  does  mean 
that  the  error  in  position  is  exactly  the  value  shovm,  rather  it  ex¬ 
presses  the  probability  that  the  true  error  in  position  will  not  be 
larger  than  the  error  given. 

Positional  error  should  be  expressed  by  precision  indexes  whicl, 
immediately  identify  the  form  and  probability  represented  by  a  given 
error.  For  example,  let  the  circular  probable  error  (CPE)  of  a  geodetic 
position  equal  100  fe“t.  Then  the  form  is  circular  •  Tile  magnitude  100 
feet  and  the  probability  (tOji  by  definition  of  CPE)  are  derived  from  a 
statistical  treatment  of  k.nown  or  estimated  error  components  comprising 
the  total  positional  error.  The  statement  Infers  a  50-50  chance  that 
the  geodetic  position  in  question  does  not  vary  more  than  100  feet  from 
the  true  geodetic  position,  v/hen  the  error  magnitude  is  increased  by  a 
statistical  factor,  greater  probability  is  achieved.  Multiplying  100 
feet  by  1.8227  yields  a  90^  probability  that  the  positional  error  will 
not  exceed  l82  feet. 

Errors  in  different  forms  are  more  easily  understood  when  precision 
indexes  common  to  linear,  circular,  and  spherical  error  distributr.ons  are 
used.  Precision  indexes  suitable  for  expressing  positional  error  in¬ 
clude  (l)  the  linear,  circular,  and  spherical  standard  errors  represent- 
ir.g  39-35°/^>  and  19.9‘l^  probabilities,  respectively,  (2)  the 

linear  probaole  error,  circular  probeble  enor,  and  spheri-’ed  probable 
error  reprosenti.ig  yjf  probability  x.i  eacn  aisunnutrcn,  (l)  the  map 
accuracy  standard,  circular  map  accuracy  standard,  and  spherical  accuracy 


56 


'Standard  ropre&entlnc  a  probauility  level,  and  (*0  a  probability 
level  ar;proachinc  near-certainty  for  each  distribution  which  the 
positional  error  is  theoretic ally  unlikely  to  exceed;  (a)  three 
si-pna  (linear,  (b)  three-i'ive  sigma  (circular, 

and  (c)  four  sigma  (spherical,  Since  error  values  are 

easily  converted  iroin  one  precision  index  to  another  in  tne  same 
distribution,  the  use  of  any  index  is  largely  a  matter  of  choice. 
Hovrever,  in  presenting  xjositional  inforn.ation,  the  positional  error 
is  best  expressed  by  either  the  90%  or  90%  probability  precision 
index  or  both. 
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^4.3.  Summary  of  Formulae  and  Conversion Factors . 


Linear  Error  Formulas 


Precision  Index  Sjcnbol 


Standard  Error 


Probable  Error 


Map  Acciu-acy 
Standard 


Near-Certainty 
Error  (Three  sigma| 


99.73^6 


’c  ''mula 


'  2 
E  (X,  -  X) 


vhere:  =  a  measured  value  of  the 

quantity  X;  X, ,  Xo  • • •  X_ 


JC  =  the  most  probable  value 
(arithmetic  mean)  of  X 
EX, 


X  3  the  error;  x  =»  -  X 

n  a  number  of  measu*'eraents 


PE  a  0.67^5  Ox 


MAS  =  1 . 64^9 


3.0000  3^ 


Linear  Error  Conversion  Factors 


50$ 

68.2756 

90$ 

99.7356 

50$ 

1.0000 

1.4826 

2.4387 

4.4475 

68.27 

0.67*45 

1.0000 

1 . 6449 

3.0000 

90 

0.*4101 

0.6&8^ 

1.0000 

1,8239 

99-73 

0.22*^ 

0  3333 

0.5483 

1.0000 

'Subscripts  denote  the  standard  error  computed  i’rom  a  sample  (o  ,  a  >  a,)- 

X  y  5; 


Circulai'  Error  Formulas 

Precision  Index 

Symbol 

Percentage 

Probability 

Formula 

Circular 

Standard  Error 

39.35^ 

ffC  =  0.5000  {Oy.  +  Oy) 

Omin/omax>  0-2' 

Circular 

Probable  Error 

CPE, 

CEP 

50^ 

CPE  =  1.177^  Oj. 

CPE  =  O.5S67  iUy  -  Oy) 
when  o^in/^max  > 

CPE-(0.2li;l  0.6621 

when  0.1  <  <  0.2  * 

CPE  -  (0.0900  0  .  +  0.6745  a  ) 
min  max ' 

when  0.0  <  <  0.]  * 

Circulai’  Map 
Accuracy  Standard 

CMAS 

CMAS  =  2.i460  Oq 

CMAS  =  1.0730  (ojj  Oy) 

when  0.2 

Circular  Near- 
Certainty  Error 
(Three-five  sigma) 

3.5d, 

99.78^6 

3.5000  0^ 

Circular  Error  Conversion  Factors 

To 

From\..^^ 

50$ 

635^ 

90^ 

99.78^ 

50 

63 

90 

99.73 

1.0000 

0,8493 

0  7071 
0.4660 
0.23?Y 

1.1774 

1 .0000 
0.8325 
0.5486 
0.3364 

1.4i42 

1.2011 

1.0000 

0.6590 

0.4040 

2.i460 

1.82?7 

1 . 5174 

l.OOOO 

0.6131 

3.5000 

2.9726 

2.4749 

1.6309 

1.0000 

'  Inhere  is  th;  minj-piara  or  smalle'f’  linsar  jtandai’d  error  of  the  two. 

‘a  circular  error  concept  is  not  ieco-jnended  for  Omin/'umax  ratios  less 
than  0.2.  However,  a  near-linear  '■)0$  probability  error  may  be  computed 
to  represent  a  CPE  for  lower  ratios  when  a  comparison  of  circular  errors 
derived  from  different  sources  is  required. 


SphericHl  Error  Formulas 


Precision  Index 

Percentage 

Probability 

Formula 

Spherical 

Standard  Error 

®s 

19.9^6 

=  0.3333(0^  “y 

W'W  >  0-35' 

Spherical 

Probable  Error 

SPE 

50^6 

SPE  =’.-5332  Os 

SPE  =  0.5127  (Ojj  <Jy  +  0^) 

when  0  .  /a  >0-35 

min'  max  ■ 

Sjherical 

Acc’uracy  Standard 

SaS 

90i 

SAS  =  2.5003  Og 

SAS  =  C.8333  (ojj  +  Oy  +  0^) 

when  a  /u  >  O.35 
min  max  — 

Spherical  Near- 
Certainty  Error 
(Four  sigma) 

99. m 

4.0000  og 

Smerical  Error  Conversion  Factors 

To 

From'''\_^ 

I9.9i> 

90$ 

6l^ 

90^ 

99. 89^ 

19.9$ 

1.000 

1.538 

1.732 

2.500 

4.000 

50 

0.650 

1.000 

1.126 

1.625 

61 

0.577 

0.888 

1.000 

1.443 

2.309 

90 

O.bOO 

C.615 

0.693 

0.433 

i.OCO 

99.89 

0.250 

0.385 

0.625 

‘a  spherical  concept 


is  not  reconunended  when  o  .  /a  it 

min  max 


less  than  0.35- 


I 
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Appendix  A 


PERCENTAGE  PROBABILITf  FOR 
STANDARD  ERROR  TNCRFJiENTS  ' 

The  following  tabic  presents  the  increments  of  linear 
circular  (on),  and  spherical  (a„)  standerd  errors  for  Intervals  of 
one  percent  probability.  Percentage  levels  corresponding  to  pre¬ 
cision  indexes  are  underlined. 

Factors  for  converting  the  error  at  one  percentage  probability 
to  another  within  the  same  distribution  are  derived  by  dividing  the 
standard  error  increment  of  the  new  percentage  probability  by  the 
standard  error  increment  of  the  given  percentage  probability.  An 
example  is  the  conversion  from  the  circular  map  accuracy  standard 
(90'r)  ■to  the  circular  probable  error  {‘^0$): 

CPE  *  1.177*^  Og 
CMAS  =  2.1‘V60  Oj. 

1.177^^ 

CPE  = - CMAS 

2.i460 

.'.CPE  =  0.5^6  CMAS 


00  0.0<.-00  o.ooor  O.OCOO 
01  0.0125  0.l4l4  0.33^39 
02  0.0251  0.2010  0.4299 


'  Reference  27. 
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37 
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39 
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Appendix  B 

THE  MOST  PROBABLE  VALUE 


Since  the  true  value  of  a  measured  quantity  is  never  known,  the 
most  probable  value  of  the  quantity  must  be  determined  from  the 
observed  values.  The  following  proof  (ref.  no.  5)  viil  •’how  that  the 
arithmetic  mean  of  the  observed  values  is  the  most  probable  value  of 
the  quantity: 


Symbols: 

X  -  an  unknown  quantity 

X,  «=■  the  observed  values  of  the  unknown  quantity; 

X^  «  Xp  Xg,  X3  ...  Xn  (1) 


the  arithmetic  mean  of  the  observed  values; 


X  - 


Ji,  X 


7  ,  or  n  X  -  Z  h 

itl  "  i»l 


(2) 


the  error  in_an  observation; 
X,  »  X.  -  X 

X  1 


(3) 


Proof: 


X.  -  X 

JL 


X  -  X 
2 


i-l  i-1 


B-1 


From  equation  (2); 

n  n  n 

t  \  •  t  h  -  t  ^  ° 

i-l  i-l  i-=l 

Tills  shows  that  the  sum  of  the  differences  about  the  mean  is  zero, 
which  was  expected,  but  if  equation  (3)  is  squared  and  then  summed: 

-  2Xj^  X  +  ("'i 

“  Xg^  -  axg  X  +  }? 


C  =  X„2.2X„X+x2 


t  -  t  ■  25c  £  *  d? 

J  _  t  *•  J  _  1  * 


The  most  probable  value  will  be  i’ound  when  ^  x.^  ■  0,or  the 

i*l 

n  P 

most  probable  value  of  X  will  be  that  which  makes  ^  "  a  minimum. 

1»1 

In  order  to  find  this  minimum,  differentiate  equation  (6)  with  respect 
to  X  and  equate  to  0: 


••  11 

^  Xj.2  =  -  2  y  X^  +  2  n  X  =  0 


•  -  AX, 

.  .  X  =  £  J:  (7) 

i»l  n 

Equation  (?)  proves  that  the  mean  value  "X  is  the  most  probable 
value  of  a  set  of  independent  observations.  Therefore,  in  the  determi¬ 
nation  of  the  residual  value  it  is  correct  to  use  the  mean  value  for 


an  approximation  of  the  true  value. 
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Appendix  C 

PFOPAGATION  OF  ERRORS 

A  quantity  Is  computed  •from  two  measured  quantities  a  and  b, 
where  ^(a/b)  denotes  a  function  of  a  and  b.  The  error  Af  of  f^  is 
affected  by  th»»  errors  In  both  a  and  b:  de.  and  /ib.  Assuming  a  and 
b  are  Independent,  and  the  errors  £»,  /5b  are  randomly  distributed,  the 
combined  error  Af  can  be  computed,  (ref.  nos.  3>  15) 

Let: 

-  f  («!,  bp 
•  f  («2.  b^) 


f  a  r  ( E  •  b  ) 

*n  '  n^  n' 

The  measured  values  of  a  and  b  may 
a  and  b.  The  most  probable  value  of  f 


(1) 

be  averaged,  obtaining  the  values 
is  f,  (from  appendix  B),  where: 


f  ■  f  (a,  b) 


ano : 


Af,  -  f. 


(2) 


In  order  to  find  the  value  of  take  the  partial  derivative  of  fj^: 


Af. 


?£i 

9a< 


Sbi 


Aa^  + 


(3) 


From  Appendix  B,  5^  Af .  -  0 

Computing  the  sum  of  the  squares  of  equation  (3); 


■  !5)  = (9  (|i  *  (9 


-.2  _  M  .  2. 


(Afn)  -  —  ^5Bn  +2  ~j 


utji  2 

+  r- 


Since : 


also: 


dag  da„ 


^  _  K 

db2  db„ 


a  constant; 


a  constant: 


Z  ■  0  Z  0 Z 

1-1  1-1  1-1 

•ig's.- 


ruvliilng:  through  by  n: 


y  fi!  .  f:i!  ts  (^1  (^y 

Z— ,j  n  i5a./  /  ■  n  *5aJ  ■Sbj^v  ■  n 

i-1  i»l  i*--! 

n 

UJ  zL  n 


By  defirition: 


V  ^  2  V’  -s?® 

L  “  %  -L 


TT  -  Cf, 


Since  a  and  b  are  independent: 


l^^i)  2_,  " 


Therefore ; 


i  2 


°a  ^  % 


Equation  (8)  is  the  general  form  for  the  propagation  of  independent 
errors,  and  can  be  expsinded  to  cover  any  number  of  quantities  (a,  b,  c,. 
d,  ....)•  It  is  imperative  that  each  element  represent  the  same  preci¬ 


sion  index  in  the  equation. 


Special  Rules  for  Error  Propagation 


Rule  1.  Sum  and  Difference;  f={a  +  b+  ...)orf=(a-b-  ...) 


Sa 


rr  »  1 

Ob 


Placing  (9)  in  the  general  equation  (8): 


(9) 


Of  =  +  0^  ?  +  .  .  .  (10) 

The  absolute  standsurd  error  of  a  quantity  computed  from  the  sum  or 
difference  of  measured  quantities  is  equal  to  the  square  root  of  the 
sum  of  the  squared  standard  errors  of  the  measured  quantities. 

This  is  the  form  most  frequent}^  encountered. 


Rule  2.  Product  of  Factors  Raised  To  Various  Powers :  f  =  a"*  b'^ 


m-1  i,q  . 

^  =  ma  b^  and  ^  ^  9^ 


Placing  (ll)  into  equation  (8): 

„  ^2m-2  r^q  „  2  ,  „2m  2  ,  2q-2  „  2 

Of  m  a  b  +  a  q  b  a^‘=- 

Dividing  through  by  f  ='\/ a^^®  b^'^; 


2  „2m-2  .  2q  ,  2 
m  a  ^  ‘^a 


a2”'  q2 


^2“  b^^ 


(11) 


(12) 


(13) 
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Hule  3-  H.mple  Product  or  Quotient:  From  the  proceeding  rule, 

f  =  a*”  let  m  =  1,  q  =  ±  1 . 

Then,  f  =  eb,  or  f  =  a/b. 

From  Equation  (13): 


where  a^/f  is  the  fractional  standard  ei-ror. 
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DERIVATIOR  AlID  SOLUTION  OF  THIS  WO-DIMENSIONAI. 
PROBABIUTy  DIOTEIBUTION  FUNCHON 

Derivation,  (ref.  no.  24)  The  probability  density  functions 
the  independent  errors  "x''  and  "y"  f',re; 


The  two-dimensional  prooability  distribution  function  is: 


2  r2  fsinS  0  oos2 


e 


2  L  0\j 


Ox 


r  dr  d© 


r  dO  dr  (small  increment  resulting 
fiom  dx  and  d:;) 


Using  identities:  sin^  ©  »  i  (l  -  cos  2©) 

2 

coa^  ©  “  -  (l  +  cos  2©) 


P(R) 


K 

I 

2*  ajc  J  J 


R  2*  2 

-  -  £1 

e 


[1  -  cos  2d  1  +  cos  29^ 


r  dr  d© 


r“0  ©“0 

Rearranging  terms: 


P(R)  - 


‘^x  ‘'y 


/ 


r"0 


R  2 


re 


r  1  1  ■ 

_  S  *  .2 
LOy  Ox  J  4 


I 

©-0 


cob2© 


Let  ^  •  20, 
dj^  -  2d©; 
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d©dr 


The:.: 


where  Iq  is  a  Bessel  Function,  z6to  order,  modified  first  Idnd. 
Therefore : 


r=0 

2 .  Special  Case  of  Two-Dimensional  Probability  Distribution  Function . 
vnien  Ojj  “  Oy  ==  Op  (ref.  nos.  l8,  24)^  from  equation  (3): 


l')-3 


Since : 


Then: 


3.  Modified  Forn;  of  the  Two-Dimenbional  Prohabiilty  Distribution 
Function,  (ref.  no.  2^)  To  solve  equation  (3)  by  the  use  o^  tables, 
the  equation  must  be  modified.  From  S.O.  Rice's  ''Properties  of  Sine 
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Wsve  Plus  N-dso"  Bell  System  Technical  Journel  V'-’l.  27  No.  1,  January, 
PP  i09”157: 


I,  (kx)  =  J 

0 


e"^  (vk)  dv 


Modifying  equation  (3): 


djj  e>y 


i.  p  p  2  * 

4a /J 

l.  Uy  Ux  J 


Step  A 
I/etting : 


:  “y 

- rr  1  +  — 7 

,  _  2  _  ' 


_lL  1  +  JL  dr 
Ucy"  I  ax"/ 


2  y 

4ay  dv  s*  2r  1  +  — ^  dr 


-  2  ^ 

r^  \ 

To  get  the  quantity  2  ■  ij  in  the  form  of  (vk): 

^a  'o  I 

t-  y  X  J 


Let  a  -  ■-  where  le  the  amaller  of  the  two: 

y 


Step  C 


Getting  Ojj  and  Oy  in  terms  of  a: 


Conblning  Steps  A,  B,  C,  D  and  equation  (3): 


0 

Revritine  equation  (9): 


X 


U-'I 


h.  Solution  of  Mcdil'ied  Function,  (ref.  nos.  'i2,  2^)  fo  compute  the 
CPE  (CPE  a  R  when  P(R)  «  0.5)  for  values  of  Oy  «  v  .b"  and  =  nTjT, 
two  methods  are  available: 

Method  1: 

To  determine  the  value  for  x  by  Rice's  table  of  Ig  (vk)  dv,  enter 
the  table  with  values  of  k  and  tne  required  probability. 

P(R)  «  509^  probability;  a  ^  =  O.C165;  a^  =  0.6667;  k  *  »  0.2 


f  Ip  (vk)  dv 

1  r  ef^J 


X, 

.50  (1  ^  f 

Oo  =  J  ^  I, 


(vk)  dv 


0.5103  * 


A 


In  (vk)  dv 


Enter  the  tables  with  k  ■  0.2  and  interpolate  for  O.5IO3  to  get  the 
value  of  X. 


.6  ^ 

^?1Y 

5103 

.8  = 

5516 

X  = 

0.71732 

X  ■ 

R^ 

ha  2 

B'1-  ”■ 


71732 
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The  radius  of  the  'yO'fo  probability  circle  (CPE)  resulting  from  Oy  is 

R  =  1.0713  Oy 


Method  2: 

Using  tables  computed  by  Arthur  Grad  and  Herbert  Solomon: 
From  equation  ( 2 ) : 


P(R)  »  P  j-v/xS  +  y2  “  *  P  I  +  y^  -  R^j 

Since  x  and  y  have  unit  standard  errors,  they  can  be  written  as: 
X  =  (Jx  ^  y  “  {jy  y* 

Therefore : 

P(R)  B  P  |oy^y^  +  Oj(® 

-  p 


y?  “5  r2 


'^y-- 


(11) 


From  Grad  and  Solomon  Tables: 

^  I®1  ^1^  ^  ^ 


P  !yo^  +  ~  -  f 

£,  Kg  ±  r.g 


(12) 


Correlation  between  equations  (ll)  and  (l2)  will  permiL  use  of  the 
tabled  values. 
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Enter  the  tables 
Then  intei’pciace 


rfith  values  of  aj,  a,^  and  ths 


foi-  values  of 


requi'"td  probability. 


Lnce 


.6 


t  «  .6  a  4559 
?000 
•7  5080 


t_ 

^2 


0.6846»t 

— JT- 


R 


2 


R  :>  ^7*6846^ 
°y  *  -6 


1.068 


R  *’  1.068  a 

y 
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Appendix  £■ 

DERIVATION  OF  THE  SPHERICAL  PROBi'^IUTy  DISTRIBUTION  FUNCTION 


Integra' injf  by  parts: 


L 


r  = 


In  order  to  use  approximation  formula  ( Mathematicai.  Tables  arid  Other 
Aids  to  Computations,  Vol.  XI,  No.  60,  October  1957,  pp  265,  "A  For;  ila 

for  the  Approximation  of  Definite  Integrals  of  the  Noniuil  Distrllution 

^  -t‘-/2 

Function"),  P{S)  mu-st  be  transformed  to  the  integral  of  e  dt. 


Lettiiig  C  ,  dS  *•  0  dC,  where  c;„  “  constant: 


■VTf- 


From  above  reference  when  x  -  0: 


2  S'’ 

(f  C=  -2-  _  C_ 

2  '  *^s  2 

Ce  +  J  e  dC 

C=0 

tjVinw  y  ^  n  ♦ 


e  dt  ~ 


X  +  0.8  e 


m 

-I 

y 

r 

-  c  e  +  1.253 


C  +  0.8 


(8) 
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